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ABSTRAGT 

Here we analyze radio, optical, and X-ray data for a peculiar cluster Abell 578. This clus¬ 
ter is not fully relaxed and consists of two merging sub-systems. The brightest cluster galaxy, 
GGPG 0719.8+6704, is a pair of interacting ellipticals with projected separation ^ lOkpc, the brighter 
of which hosts the radio source 4G +67.13. The Eanaroff-Riley type-H radio morphology of 4G +67.13 
is unusual for central radio galaxies in local Abell clusters. Our new optical spectroscopy revealed 
that both nuclei of the GGPG 0719.8+6704 pair are active, albeit at low accretion rates corresponding 
to the Eddington ratio ^10“^ (for the estimated black hole masses of ^ 3 x 10^ Mq and ^10^ ^©)- 
The gathered X-ray {Chandra) data allowed us to confirm and to quantify robustly the previously 
noted elongation of the gaseous atmosphere in the dominant sub-cluster, as well as a large spatial 
offset 60kpc projected) between the position of the brightest cluster galaxy and the cluster cen¬ 
ter inferred from the modeling of the X-ray surface brightness distribution. Detailed analysis of the 
brightness profiles and temperature revealed also that the cluster gas in the vicinity of 4G+67.13 is 
compressed (by a factor of about ~ 1.4) and heated (from 2.0keV up to 2.7keV), consistent with 
the presence of a weak shock (Mach number ^ 1.3) driven by the expanding jet cocoon. This would 
then require the jet kinetic power of the order of ^ 10^^ ergs“^, implying either a very high efficiency 
of the jet production for the current accretion rate, or a highly modulated jet/accretion activity in 
the system. 

Subject headings: galaxies: active — galaxies: individual (4G +67.13 ) — galaxies: clusters: individual 
(Abell 578) — intergalactic medium — galaxies: jets — X-rays: galaxies: clusters 


1. INTRODUCTION 

Relativistic jets produced in active galactic nuclei 
(AGN) are thought to play a major role in shaping the 
co-evolution of supermassive black holes (SMBHs) and 
galaxies via widely discussed, though still hardly un¬ 
derstood, feedback processes (see the topical reviews by 
Eabian 2012; Kormendy & Ho 2013). The most direct 
manifestation of the impact the jets make on their en¬ 
vironment is the presence of cavities inflated in the in- 
tergalactic/intracluster medium at the position of jet co¬ 
coons (or ‘radio lobes’), which are often accompanied by 
shock waves driven in the surrounding hot. X-ray emit¬ 
ting gas (see, e.g., McNamara & Nulsen 2012, and refer¬ 
ences therein). 

The more luminous the jets are, the more dramatic 
their influence is on the ambient medium. It has been 
argued in particular, that high-power jets are necessary 
to explain the observed properties of the central parts 
of rich clusters as due to mechanical gas heating by ex¬ 
panding lobes (see, e.g., Vo it & Donahue 2005; Math¬ 
ews & Guo 2011). These high-powers are equivalent 
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^ Based on service observations made with the WHT operated 
on the island of La Palma by the Isaac Newton Group in the 
Spanish Observatorio del Roque de los Muchachos of the Insti¬ 
tute de Astrofisica de Canarias. 


to those in radio-loud quasars and their parent popu¬ 
lation of Eanar off-Riley type-H radio galaxies (hereafter 
ERHs). Yet, most of the jetted AGN found in the centers 
of nearby clusters are low-power systems, with large-scale 
edge-dimmed ERI radio morphologies (e.g., Zirbel 1997; 
Wing & Blanton 2011). One notable exception from this 
rule is the ERH radio galaxy Gygnus A, located in the 
center of a rich cluster (Garilli & Barthel 1996). This 
source is considered to be representative of high-redshift 
radio galaxies since in the earlier Universe, when the bulk 
of the feedback process took place, powerful AGN with 
luminous jets are being found also in dense environments 
(e.g., Siemiginowska et al. 2005; Belsole et al. 2007; An- 
tognini et al. 2012) 

In order to increase the sample of ER H radio sources 
studied in nearby clusters, we examined the NRAO^ Very 
Large Array (VLA) imaging survey of Abell clusters with 
richness R>0 (Owen et al. 1992; Owen & Ledlow 1997). 
We found that within the redshift range 2 : < 0.25 and 
at the northern declinations visible to the VLA, only 
two ERH radio galaxies are hosted by brightest clus¬ 
ter galaxies (BGGs) in this sample: PKSB1358—113 in 
Abell 1836 and 4G+67.13 in Abell 578. Eor these two 

^ The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under cooperative agree¬ 
ment by Associated Universities, Inc. 


2 



Fig. 1.— VLA 4.7 GHz contours (levels start at 0.2 mjy beam“^ and increase by factors of two) of the two radio galaxies found in the 
Abell 578 cluster superimposed on smoothed optical SDSS z-band images. The left panel shows the BCG-hosted 4C+67.13, and the right 
panel shows the WAT radio galaxy located in the southern part of the cluster. The fainter second galaxy adjacent to the 4C+67.13 host, 
which is seen through the eastern radio lobe (in projection), is the other member of the galaxy pair CGPG 0719.8+6704. A bright optical 
source visible to the south of 4C+67.13 is a foreground star. Throughout, coordinates are in J2000.0 equinox. 


systems, we have re-analyzed archival VLA maps, ac¬ 
quired new optical spectra using the William Herschel 
Telescope (WHT), and obtained new X-ray data with 
Chandra and XMM-Newton. In this paper we present 
the analysis^ of 4C +67.13/Abell 578, which is a peculiar 
example of a cluster in formation. The companion study 
of PKSB1358—113/Abelll836 is presented in Stawarz et 
al. (2014). 

Based on a detailed analysis of the velocity mea¬ 
surements for the cluster member galaxies, along with 
ROSAT PSPC data, Gomez et al. (1997) concluded that 
Abell 578 (cluster richness R = 0) is not relaxed, and 
consists of two interacting sub-clusters with a velocity 
difference of ^ 1300 km s“^. The BCG is associated with 
the dominant cluster, and is the brighter of a pair of in¬ 
terconnected red spherical compacts as noted by Zwicky 
(1971) and catalogued by him as CGPG 0719.8+6704. 
The radio galaxy 4G +67.13 is hosted by the BGG and is 
characterized by an FRII large-scale morphology (Owen 
& Ledlow 1997, designated as 0719+670 therein), and 
hence potentially by a high jet kinetic power. At the 
southern end of this dominant cluster is the sub-cluster 
containing another radio galaxy, designated as 0720+670 
by Owen & Ledlow (1997), which is of the ‘wide-angle 
tailed’ (WAT) type (see Section 2). Both radio galaxies 
are at rest with respect to their sub-clusters (Gomez et 
al. 1997). 

Below, we assume a standard cosmology with Hq = 
71 kms“^ Mpc“^, = 0.27 and Oa = 0.73, so that 
the redshift of the target 2 ; = 0.0866 (Owen et al. 1995) 
corresponds to the luminosity distance of dp = 391 Mpc 
and the conversion scale of 1.604 kpc per arcsecond. 

2. RADIO DATA 

Archival VLA data at two frequencies were analyzed to 
compare the 4C+67.13 radio morphology at arcsecond- 

^ The XMM-Newton data we obtained in this case was entirely 
contaminated by background flares and is not considered in this 
paper. 


resolution with the newly acquired Chandra X-ray maps. 
We combined multi-configuration 1.5 GHz data obtained 
in 1991 when the VLA was in its A- (Sep 8; program 
AL238, 6 min. integration) and G-arrays (Jan 29; AO104, 

5.5 min.), together with B-array data from 2001 May 
7 (AC582, 1.5 hrs). At 4.7GHz, we analyzed a single 
deep (2 hr) C-array dataset from 2001 July 20 (program 
AC582). We used AIPS for the standard calibration and 
DIFMAP (Shepherd et al. 1994) for self-calibration and 
imaging. The resultant map beamsizes were 3.12" x 2.03" 
(position angle, PA = 43.2°) at 1.5 GHz and 5.05" x 
3.21" {PA = 12.4°) at 4.7 GHz. 

The VLA contours of 4G +67.13 superimposed on the 
Sloan Digital Sky Survey (SDSS; York et al. 2000) op¬ 
tical image are shown in Figure 1 (left panel). Note 
that despite its edge-brightened FRII radio morphol¬ 
ogy, the studied radio galaxy does not obviously show 
compact high surface brightness hotspots at the edges 
of the radio lobes. With a total linear size of 47" ^ 
75 kpc and axial ratio 2, which is typical for FRHs, 
the entire volume occupied by the radio lobes is V — 

1.6 X 10^^ cm. The radio core of 4G+67.13 coincides 
with the BCG, which is the brighter member of the 
galaxy pair CGPG 0719.8+6704, and which is offset by 

38" 61 kpc from the position of the cluster center (as 

inferred from the modeling of the cluster X-ray surface 
brightness distribution; see §4.1.1 below). 

The radio core of 4C +67.13 is distinguished by its in¬ 
verted radio spectrum, with 1.5 GHz peak of 8.4 mJy 
beam“^ and 4.7GHz peak of 12.1 mJy beam“^. Sub¬ 
tracting the estimated contribution of the core flux from 
the lobes, the 1.5 GHz fluxes are 0.314 Jy (eastern) and 
0.331 Jy (western), and the 4.7GHz fluxes are 0.120 Jy 
(eastern) and 0.127Jy (western). The resultant spec¬ 
tral indices for both lobes are 0.84 ± 0.12, assuming 
10% uncertainties in each flux measurement, and the 
monochromatic lobes’ radio power reads as Li.sghz — 
1.8 X lO^^ergs”^. The total fluxes and the lobe spec¬ 
tra from the VLA maps compare well with those derived 
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Fig. 2.— The rest frame optical spectrum of 4C +67.13, obtained 
with WHT. The observed spectrum is plotted in black while the 
STARLIGHT model is drawn in red. Unremoved telluric lines are 
marked with 0. 

from integrated measurements at 1.4 GHz (NVSS Con¬ 
don et al. 1998) and 4.85 GHz (Gregory & Condon 1991) 
as well as the spectral index a = 0.88 ± 0.08 measured 
using 2.7, 5, and 10.7 GHz Effelsberg data (Reich et al. 
2000). At lower frequencies, we derived a = 0.64 ± 0.04 
by combining integrated fluxes at 74 MHz (Lane et al. 
2012), 178MHz (Gower et al. 1967), and 365 MHz (Dou¬ 
glas et al. 1996, assuming 15% uncertainty), together 
with the NVSS measurement, indicating flatter spectra 
below ^ 1 GHz. 

The right panel in Figure 1 presents the VLA contours 
of the WAT radio galaxy located in the southern part 
of the Abell 578 cluster. The distance between the radio 
core of this source and the position of the cluster X-ray 
centroid is about 158", corresponding to 253 kpc. 

Although 4C +67.13 lacks clearly defined hotspots, and 
its radio luminosity is intermediate between FRI and H 
sources according to the often considered luminosity di¬ 
vide, it should nonetheless be classified as an FRH ra¬ 
dio galaxy. The original classification by Fanaroff & Ri¬ 
ley (1974) is based on the ratio of the distance between 
the regions of the highest brightness on opposite sides 
of the radio core, to the total extent of the radio struc¬ 
ture; sources with this ratio greater than 0.5 are classified 
as FRHs (“edge-brightened radio morphology”). In the 
case of 4G +67.13, we obtain the ratio of 0.65 for the 
higher resolution map from Owen & Ledlow (1997), and 
0.58 for the lower resolution map (Figure 1, left). The rel¬ 
atively low radio luminosity of the analyzed object com¬ 
pared with ‘classical’ FRHs should not be considered as 
a surprise, as recent studies reveal the existence of a large 
number of FR Hs near and below the classical luminosity 
divide (Best 2009; Gendre et al. 2013). An alternative 
interpretation could be that the edge-brightened appear¬ 
ance of the radio lobes in 4C+67.13 is in fact only due 
to a special alignment of a ‘narrow-angle tailed’, i.e. in¬ 
trinsically edge-dimmed radio structure viewed with its 
prolonged tails along our line of sight. That 4C +67.13 
was selected for study from a parent sample (Owen & 
Ledlow 1997) where the BGGs predominantly show clear 
FR I morphology, would indicate an increase in the prob¬ 
ability for such a special alignment. Although we think 
it is improbable that both tails would be aligned in such 
a way that their extended diffuse radio emission is un- 
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Fig. 3.— Observed [Nil] A6548,6584 + Ha blended emission lines 
of 4C +67.13 (plotted in black), with Gaussian profiles fitted to the 
emission lines overlaid in red. The reconstructed blend is drawn in 
blue in the top panel while the residuals are shown in the bottom 
panel. 

detectable, and the resulting structure — characterized 
notably by sharp boundaries — appears highly symmet¬ 
ric in projection, we can not completely rule out this 
scenario. 

3. OPTIGAL SPEGTROSGOPY 

Optical spectroscopic observations of the BCG host¬ 
ing 4C +67.13 and its neighboring galaxy were obtained 
with the WHT on 2010 April 21. Longslit spectra were 
taken using the Intermediate dispersion Spectrograph 
and Imaging System (ISIS) double-armed spectrograph 
with a chosen slit width of 1" (corresponding to 1.6 kpc). 
The slit was centered on the nucleus of 4C+67.13 and 
was oriented along the position angle of the galaxy pair 
{PA = 93°) so data were obtained for both objects. To 
identify and reject cosmic rays, the integration time was 
split into two 900 s exposures per arm. The average see¬ 
ing FWHM during the observing run was about 0.8". 

In the blue arm, the grism R300B was used and the 
detector was a thinned, blue-sensitive EEV12, array of 
4096 X 2048 (13.5 micron) pixels with a spatial scale of 
0.20"pixel“^. The binning factor was 1x1 which yielded 
a wavelength coverage between about 3300 A and 5300 A 
with a dispersion of 0.86 ApixeR^. The instrumental 
resolution was 3 A (FWHM) corresponding to ainst ^ 
90kms“^. In the red arm, the grism R158R with an 
order blocking filter GG495 was used. The detector was 
the default chip for the ISIS red arm, RED+, a red- 
sensitive array of 4096 x 2048 (15.0 micron) pixels with 
spatial scale 0.224" pixeR^. The binning factor was 1 x 
1 which yielded a wavelength coverage between about 
5300 A and 9000 A with a dispersion of 1.82ApixeR^. 
The instrumental resolution was 6 A corresponding to 
c^inst ^ 102kms“^. 

Reduction steps were performed for both spectral 
ranges separately using the NO AO IRAF packages. A 
master bias frame was created by averaging all the bias 
frames obtained during the observing night and sub¬ 
tracted from the science frames. Also, a master flat 
field frame was created and all 2D science frames were 
corrected for flat field, then cosmic rays were removed 
in combining the science exposures. Accurate wave¬ 
length calibration was performed using ArNe lamp ex¬ 
posures obtained before and after every target exposure 
and checked using sky lines, and the correction for op¬ 
tical distortion was applied. The contribution from the 
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TABLE 2 

Results of the 2-dimensional fitting to the X-ray image 
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Fig. 4.— The exposure-corrected Chandra 0.5-7.0keV image of 
Abell 578 with the VLA 4.7 GHz contours superimposed. X-ray 
point sources, including the nuclei of the two cluster radio galaxies, 
are removed from the image. The pixel size is equal to 4 times of 
the ACIS pixel size (i.e., 4 x 0.492" = 1.968") and the image is 
smoothed by a gaussian with cr = 16 X the ACIS pixel size. Open 
black triangle and square denote the position of the cluster center 
obtained from the 2-dimensional fitting with elliptical and isotropic 
2-D /3 models, respectively (see §4.1.1.) 

sky was determined from the sky regions at the two sides 
of the resulting spectrum, and then subtracted. The ex¬ 
traction of the ID spectra was performed using APALL 
task. Scientific exposures were flux calibrated using the 
standard SP1257+038. The slit position angle falls along 
the radio jet/lobe axis, which is perpendicular to the op¬ 
tical photometric major axis so light from a region of 
14.3" X 1" (23kpcx 1.6 kpc) was summed into the ID 
spectrum. 

Before fitting the stellar continuum, the ID spectra 
were corrected for Galactic extinction using the Schlegel 
et al. (1998) maps and the extinction law of Cardelli 
et al. (1989), shifted to the rest frame using a redshift 
of z = 0.087227 (Gomez et al. 1997) and resampled to 
AA = 1 A. Analysis of the galaxy stellar continuum was 
performed using the STARLIGHT code (Gid Fernandes 
et al. 2005) which fits an observed spectrum with a lin¬ 
ear combination of the simple stellar populations (SSPs) 
extracted from the models of Bruzual & Gharlot (2003). 
Bad pixels and emission lines were masked and left out 
of the fits. Two bases of SSPs used in these fits: (1) 
= 150 with 25 ages between 1 Myr and 18 Gyr and 
six metallicities from Z = 0.005^0 to 2.5^0, and (2) 
= 45 with 15 ages between 1 Myr and 13 Gyr and 

TABLE 1 

Optical emission lines in the 4C+67.13 nucleus 


Line 

Flux 

[10“^^ ergcm“^ s“^] 

Position^ 

[A] 

FWHM 

[kms“^] 

[NII]A6548 

0.30 

6545.6 

485.7 

Hoinarrow 

0.32 

6561.7 

484.7 

Hcq^road 

2.12 

6567.9 

2336.9 

[NII]A6584 

0.895 

6581.0 

485.7 


^The redshift of 2 ; = 0.087227 (Gomez et al. 1997) is adopted. 


Parameter Value 



Gluster region 
isotropic 

elliptical 

Gore radius rc [kpc] 
Genter position 
Genter position 5^ 
Ellipticity e 

Ellipticity 6 [degree] 
Amplitude^ [10“^] 
Index ^ 

reduced G-stat/dof 

118+m 

7''25™0P.4l°;® 

+66°59'10'T3 

Q oq-|-0.30 
Z.OC)_o 25 

U.05_o 10 

1.17/600 

ITOTI 

7*'25“00=.6lo;g 

+66°59'09"l2 

U.OZ_o 05 
q 9+4.3 
o.^_4_2 

^.au_o 26 

u.Dy_o 12 
1.11/598 

Background region 

Gonstant^ [10 ^] 

1 54+U.U2 

-‘-•^^- 0.02 

reduced G-stat/dof 

1.15/1950 


^12000 equinox 
^photon pixel“^ cm“^ s“^ 


three metallicities from 0.004 to 2.5^0. The results of 
the stellar populations and velocity dispersions obtained 
from the fits made with these two bases were consistent. 

Stellar population fits to the observed spectrum of the 
BGG (Figure 2) revealed only old populations with 11 
and 13 Gyr and metallicity of O.O2Z0; no traces of young 
populations (with ages below a few Gyrs), indicative of a 
recently enhanced starformation, have been found in the 
modelling. The velocity dispersion obtained from the 
fits was cr* = 321 ± 17kms“^ which gives a SMBH mass 
of logMBH/Al0 = 8.96 ± 0.22 based on the Mbh — 
relation from Tremaine et al. (2002). Taking instead the 
Mbh — relation from Ferrarese & Ford (2005), one 
obtains logMBH/Al0 = 9.21 ± 0.27. 

The [Nil]A6548,6584 and Ra lines in the stellar light 
subtracted spectra were fitted simultaneously with Gaus¬ 
sian profiles. Fitting was performed using the SPEGFIT 
task in the STSDAS external IRAF package. Lines of the 
same ion were assumed to have the same offset and width 
and the additional constraint [Nil] A6584/[NII]A6548 = 3 
was further imposed on the line flux ratio. The results 
are presented in Figure3 and Tablet (where the pro¬ 
vided line intensities are not extinction corrected). De¬ 
tailed optical activity classification using line diagnostic 
diagrams (Kewley et al. 2006; Buttiglione et al. 2010) 
is not possible in the case of 4C+67.13 because only 
[Nil] and Ra emission lines were detected with S/N ra¬ 
tio larger than 3. Nevertheless, these lines can still be 
used for a spectral diagnostic following Gid Fernandes 
et al. (2011), indicating that the 4G+67.13 nucleus is 
of a LINER type (‘low-ionization nuclear emission-line 
region’). The Ra luminosity I/Ha = 5.9 x 10^^ergs ^ 
translates to the bolometric accretion-related luminosity 
Lnuc — 10^^ ergs“^ (Sikora et al. 2013), or in the Edding¬ 
ton units A = I/nuc/AEdd — 10“"^. This low accretion 
rate is confirmed by the analysis of the Chandra data for 
the 4C+67.13 nucleus presented below in §4.2.2. 

The analogous analysis performed for the second, 
fainter galaxy of the GGPG 0719.8+6704 pair revealed 
= 233 ± 97kms“^, log MBH/Af0 8.3 — 8.7, old stel¬ 
lar populations with a small addition of a younger one 
(1.43 Gyr, Z = 0.05 Zq)^ and a LINER-type nucleus with 
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Fig. 5. — Chandra zoom-in images of the central part of the Abell 578 cluster around the position of 4C +67.13 (left), and of the southern 
part of the cluster around the position of the WAT radio galaxy (right), with the VLA 4.7 GHz radio contours superimposed. X-ray point 
sources are not removed from these figures. The pixel size is equal to a half of the ACIS pixel size (i.e., 0.5 x 0.492''). The image is 
smoothed by a Gaussian with a = AGIS pixel size. 


Lhcx — 10^^ergs ^ or A 5 x 10 

4. GHANDRA X-RAY OBSERVATIONS 

Abell 578 was observed with Chandra on 2010 May 
29 (obsid=11749) and July 23 (obsid=12225) using the 
ACIS-S detector. The source was placed on the back- 
illuminated CCD (S3) and the observations were made 
in very faint (VFAINT) mode with a total exposure time 
of 39.3 ks. 

Analysis was performed with the Cl AO version 4.5 
software. We processed the data by running Cl AO 
tool chandra_repro and applied the calibration files from 
CALDB 4.5.5.1. In order to obtain the highest angular 
resolution image, the sub-pixel algorithm EDSER was 
also applied. Event files of the two observations were 
merged into a single event file with reproject_obs for the 
image analysis. Point sources detected with wavdelect 
were removed, including also the nucleus of 4C+67.13, 
whose spectrum is analyzed separately in §4.2.2. Be¬ 
sides the 4C+67.13 nuclear source, there were 15 X-ray 
sources found in the cluster region (< 2.5' of the cluster 
center obtained in §4.1.1) and 30 sources in the ana¬ 
lyzed region. All events corresponding to the detected 
point sources are excluded from the analysis of the ex¬ 
tended emission. Spectral extraction was performed with 
specextract and modeling was done in Sherpa (Ereeman 
et al. 2001; Refsdal et al. 2009). The spectra and re¬ 
sponse files were extracted from the individual data and 
merged with combine_spectra. All spectral models were 
fitted to the data in the 0.5-7.0keV energy range. Except 
for the radial profile analysis, we used the C-stat fitting 
statistics (Cash 1979) and the Nelder-Mead optimization 
method (Nelder & Mead 1965). 

4.1. Image Analysis 

Eigure4 shows the exposure-corrected and smoothed 
Chandra image of Abell 578 in the 0.5-7.0keV energy 
band. Events of the two separate observations were 
summed to create the image. The green contours are 
from the VLA 4.7 GHz radio map discussed in § 2 above. 

Eigure5 presents the Chandra zoom-in images of the 


central and southern parts of the cluster, around the po¬ 
sitions of the two cluster radio galaxies (left and right 
panels, respectively). Point sources were not removed 
from these figures. No obvious small-scale X-ray sub¬ 
structures adjacent to the radio lobes can be noted, most 
likely due to the limited photon statistics. 

4.1.1. 2-Dimensional Fitting 

To estimate the position of the cluster center and to 
extract the general characteristics of the cluster large- 
scale morphology, two-dimensional fits to the Chandra 
data was performed. Eigure6 shows the corresponding 
images of the background regions (upper panel) and of 
the cluster (middle and lower panels). Total counts are 
11720 for the background and 8058 for the source. 

Eirst, we modeled the background image with the con¬ 
stant model (const2d) and the source image with the 
isotropic 2-D [3 model (beta2d). The regions around 
the two radio sources were removed from the analy¬ 
sis; these correspond to the circular black spots in the 
cluster images in Eigure6. In the fitting, an expo¬ 
sure map created by running flux_obs was used to cor¬ 
rect non-uniformity of the effective area in the ACIS 
chip. The background model was fixed when fitting the 
source region. The fitting results are listed in Table 2. 
In particular, we found the position of the cluster cen¬ 
ter ((a,J) = (7^25^01®.dlg’g,+66°59'10"^3), the cluster 
core radius rc = 118^29 index (3 = 0.581q‘{q 

(reduced fit C-statistic 1.17/600), consistent with the 
values obtained by Gomez et al. (1997). 

We also modeled the source image with the elliptical 2- 
D p model. The results are given in the third column of 
Table 2. The thus evaluated position of the cluster cen¬ 
ter is consistent with the result of the isotropic (3 model, 
and the ellipticity parameters are e = 0.321 q;q 5 and 
0 = —3.2l4 2deg, while the other parameters changed 

only slightly {rc = 170l3gkpc, [3 = 0.69lo;i2)- The 
elliptical model gives, however, a better fit with the re¬ 
duced fit C-statistic of 1.11/598. In this way, we con¬ 
firmed and robustly quantified the E-W elongation of 
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Fig. 6. — Two-dimensional analysis of the background (upper panel) and the Abell 578 cluster region (middle and lower panels). The 
ellipticity parameters are fixed to the circular symmetric values in the middle panel, while they are free parameters in the lower panel. 
The images of the data (left), adopted beta2d+const2d model (center), and residual (right) are shown in the same color scale. The green 
contours overlaid on the residual images represent ±2-sigma of the error of the data count in each pixel. The black regions in the images 
are ignored in the fitting. The pixel size is 9.84" (20 times of the original ACIS pixel size). 


the cluster mentioned by Gomez et al. (1997). 

The positions of the cluster center obtained from our 
analysis are shown in Figured where the black triangle 
and square denote the values obtained with the elliptical 
and isotropic models, respectively. The cluster center po¬ 
sition is clearly offset from the position of the 4C +67.13 
core. Our radio data indicate that the position of the 
4C +67.13 core is (7^24^54M50, +66°59T0.30"). There¬ 
fore, the offset in projection is Vx — 38", corresponding 
to 61 kpc for the elliptical cluster model. 

4.1.2. Radial Profiles 

A radial surface brightness profile of the central parts 
of the Abell 578 cluster was constructed from a pie re¬ 
gion enclosed by the magenta solid line shown in Fig¬ 
ure 7. The center of the pie region was set at the position 
of the cluster center obtained from the two-dimensional 
fitting discussed in the previous section. We extracted 
surface brightness from 20 annuli covering a radius of 
150" 241 kpc) by utilizing the dmextract tool. A pie 

region embedding within the 4C +67.13 radio source was 
ignored in order to exclude the expected contribution 
from the cluster gas interacting with the expanding radio 
lobes, as well as the X-ray emission of the radio galaxy it¬ 
self (including the non-thermal emission of the extended 
lobes as well as the radiative output of the AGN and 
circumnuclear gas). We subtracted the surface bright¬ 


ness of the region defined by a polygon as a background. 
Hereafter we refer to this background region as “back¬ 
ground A.” We fitted the background-subtracted surface 
fiux profile with the betald model using the chi2datavar 
statistics. The fitting results are presented as the ma¬ 
genta curve in Figure 8, and the best-fit parameters are 
listed in Table 3. 

In order to investigate the impact of the central radio 
source on the surrounding cluster gas and to further ex¬ 
amine the cluster elongation, we also constructed radial 
surface brightness profiles separately for the north, east, 
south, and west parts of the cluster center. The respec¬ 
tive black, red, green and blue dashed regions in Figure 7 
were used to construct these profiles. Each region was di¬ 
vided into 8 annuli, with the radio lobes and X-ray point 

TABLE 3 

Results of the X-ray radial profile analysis 


Parameter Value 

Core radius rc [kpc] 

Amplitude^ [10“^] 2.80 ^q'2q 

Index /3 0.71 ^q‘^® 

xi/doi _1-25/17 


^photon pixel ^ cm ^ s ^ 
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30.0 7:25:00.0 30.0 

Right ascension (J2000) 


24:00.0 


Fig. 7. — Chandra X-ray image of Abell 578 (gray scale) over¬ 
laid with the VLA 4.7GHz contours (yellow). Different regions 
selected for constructing radial surface brightness profiles of the 
cluster are denoted by magenta, red, black, green, and blue con¬ 
tours (see §4.1.2). Region defined by a white polygon was used to 
extract the background surface brightness (“background A”). 



Fig. 8.— One-dimensional radial X-ray surface flux profiles of the 
Abell 578 cluster extracted from different regions defined in §4.1.2 
(color-coded symbols denote the corresponding regions in Figure 7). 
Thick magenta curve represents the best-fit betald+constld model 
applied to the region enclosed by the magenta solid line in Figure 7. 


sources removed. The profiles cover the regions from the 
center to a radius of 144" 231 kpc) within an angle of 

90°. The results are shown in Figure8. The cluster sur¬ 
face brightness in the western region (blue filled circles), 
where the central radio source 4C +67.13 is located, dis¬ 
plays a clear excess over the average model. This may 
be due to the gas compression/heating by the expanding 
radio lobes. Also, for radii > 80", the surface brightness 
in the northern region (black crosses) is slightly lower 
than the average model, while that in the eastern region 
(red diamonds) slightly higher, consistent with the E-W 
elongation of the cluster atmosphere. 

The last two data points in Figure 8 for the southern 
region (green triangles) overlap with the position of the 


WAT radio galaxy. While we expect the elongation of 
the cluster to reduce the surface brightness along this 
direction compared to the average model, this effect is 
compensated for by the cluster substructure discussed 
also in Gomez et al. (1997), therefore no surface bright¬ 
ness decrement is seen. 

4.2. Spectral Analysis 
4.2.1. Cluster Gas Around 4C-1-67.13 

The analysis presented in the previous section provided 
some indication for the interaction of the 4C +67.13 with 
the surrounding gas. In order to investigate in more de¬ 
tail the impact of the radio source on the cluster environ¬ 
ment, we extracted spectra from four smaller pie regions 
around the cluster center, with the outer radius of 75", 
as shown in Figure 9. We refer to the western, eastern, 
northern, and southern regions as “source 1”, “source 2”, 
“source 3”, and “source 4”, respectively. From source 1, 
the emission from the positions of the radio lobes and 
the core of 4C +67.13 were once again removed. For the 
background, we chose the polygon region A, the same as 
the one used in constructing the cluster radial profiles 
(see § 4.1.2), and a large annular region “background B.” 
While background A represents a uniform background 
in the ACIS chip — i.e., consisting of the Cosmic X-ray 
Background (CXB), Galactic Halo (GH) emission, and 
the instrumental background — background B was cho¬ 
sen to characterize the emission from the outer parts of 
the cluster adjacent to the selected source regions 1-4. 

First we fitted the spectra of backgrounds A and B 
(total counts in the 0.5-7.0keV energy band 3464 and 
3487, respectively). For the instrument background, we 
used the model defined in “acis_bkg_model.py”^ that con¬ 
sists of two power-laws and six gaussians. Parameters 
in the model were fixed at the values for ccd_id = 7, 
and only the normalization was a free parameter. We 
used a constant ARF by setting all values in specresp 
column to 100. The CXB and GH components were 
modeled by an absorbed power-law model and an unab¬ 
sorbed thermal plasma model APEC, respectively. The 
photon index of the CXB was fixed at 1.41 (De Luca & 
Molendi 2004) and the temperature of the GH was fixed 
at 0.2 keV (Kuntz & Snowden 2000). The absorption col¬ 
umn density for the CXB was fixed at the Galactic value, 
4.39 X 10^^ cm The spectrum of background B was 
modeled by the background A scaled for the extraction 
region and an APEC model with Galactic absorption. 
The redshift of the APEC model for background B was 
fixed at 0.0866, and the abundance was allowed to vary. 
The obtained best-fit parameters of backgrounds A and 
B are listed in Tabled. The resulting normalization of 
the CXB component 12.4 ± 1.5phkeV“^ cm“^ s“^ str“^ 
at 1 keV is consistent with that of De Luca & Molendi 
(2004). These background model parameters were then 
fixed in the following analysis. 

In the above, we have assumed that background A does 
not contain any significant cluster emission. We verified 
this assumption by adding an APEC component with 
Galactic absorption to the model for background A and 
obtained an upper limit for the APEC normalization, 

^ https://github.com/taldcroft/datastack/blob/master/ 

acis_bkg_model.py 












Fig. 9. — Chandra X-ray image of Abell 578 (gray scale) overlaid 
with the VLA 4.7 GHz contours (yellow). Different regions selected 
for the cluster spectral analysis are denoted by red, black, green, 
and blue contours (see §4.2.1). Regions defined by white polygons 
were used to extract the background spectra (“background A” and 
“background B”). 


with the 90% confidence intervals of 7.6 x 10“^ for 2keV 
and 2.4 x 10“^ for 1 keV, where the temperature of the 
APEC component is fixed. Since the backscale value of 
source 1 is ^ 0.1, this upper limit corresponds there¬ 
fore to a negligible amount of the APEC component for 
source 1 spectrum. 

Next we used the Cl AO Sherpa extension package de- 
project^ (Siemiginowska et al. 2010) for spectral anal¬ 
ysis of the selected source regions. Here, we assumed 
spherical symmetry for the cluster gas. In the anal¬ 
ysis, the spectral model was composed of a linear 
weighted combination of a source model, the background 
B model, and the background A model was scaled by the 
backscale value: rsp(Wsrc*src_model+Wbkg*bkgB_model) 
+ bkg_scale*(bkgA_rsp(bkgA_model)). The weight was 
defined as a volume of extraction region normalized 
by a volume of a sphere enclosing background B. The 
weight was calculated in the deproject tool except for the 
source 1, for which the weight was calculated by Monte- 
Carlo integration. The weight of the source 1 region is 
0.0245852 ± 0.0000008 and that of background B for the 
source 1 analysis is 0.0462102 ± 0.0000009. We assumed 
an APEC model with Galactic absorption for each spec¬ 
trum. We fixed the abundance at 0.0, 0.5 and 1.0 at first, 
and then set it free. 

X-ray spectra and spectral fits for the source 1 and 
the combined source 2+3+4 regions are shown in Eig- 
ure 10; the model parameters for the abundance A = 0.5 
are listed in Table 5. The total counts in the 0.5-7.0keV 
energy band were 838 for soured, 884 for source2, 880 
for source 3, 773 for source4, and 2537 for source 2+3+4. 
The corresponding contour plots of temperature versus 
normalization are shown in Eigure 11. Eor all the consid¬ 
ered abundance values, source 1 (blue) is located at the 
upper right in the contour plots, which implies that the 
cluster gas around 4C +67.13 is indeed denser and hotter 

^ http://cxc.harvard.edu/contrib/deproject/ 


TABLE 4 

Best-fit model parameters for backgrounds A and B 


Model component 


Parameter 

Value 

background A 

Instrumental background 

normalization 

31.42+^;^^ 

Cosmic X-ray background 

photon index 
norm^ [ 10 - 6 ] 

1.41 (fix) 

Galactic Halo 


kT [keV] 
norm^ [ 10 “®] 

0.2 (fix) 

[T 0 + 3.2 
^•'^-3.1 

Fit statistics 


reduced G-stat/dof 

0.64/60 

background B 

Cluster gas emission 


Abundance 

0 34+u-^3 



kT [keV] 
norm^ [ 10 “^] 
density [ 10 “^ cm“^] 

Z.OU_o 37 

'-'•^'^-0.07 

n 70+0-06 

'®-0.05 

Fit statistics 


reduced G-stat/dof 

1.04/51 


^Normalization of the power-law model in units of 
photon cm“^ at 1 keV. 

^Normalization of the APEC model defined as 
10-i4/(47r [Da{1 + z)f) f neUHdV. 

than that in the other selected regions, consistent with 
shock heating due to the expanding radio lobes. The re¬ 
duced C-stat of the spectral fit for source 1 were 1.44, 
1.33, 1.36 and 1.37 for the abundance 0, 0.5, 1 and free, 
respectively. Also, for the combined source 2+3+4, the 
reduced C-stat is best at the abundance of 0.5; hereafter 
we adopt the fitting results with the abundance fixed 
at this value. The temperature ratio of source 1 to the 
combined source 2+3+4 is then 1.37lo.i85 
sity ratio is 1.36lo!o6* These values correspond to shock 
Mach numbers of 1.38 and 1.24, respectively. 

4.2.2. 4C+67.13 Nucleus 

Einally, we analyzed the X-ray spectrum of the active 
nucleus in 4C +67.13. The spectrum was extracted from 
a circular region with 1.5" radius while the background 
spectrum was extracted from the annular region from 3" 
to 11". The total number of counts in the 0.5-7.0keV 
energy band is 59 for the source region and 88 for the 
background region. 

The background spectrum was fitted by the model used 
in the fitting of background A plus an APEC model mod¬ 
erated by the Galactic absorption. The normalizations 
of the instrumental background, as well as of the CXB 
and GH, were fixed at the values obtained in section 
4.2.1, and scaled by the areas of the extraction regions. 
Abundance was fixed at 1.0, and the redshift was fixed 
at 0.0866. The best-fit model parameters obtained are 
listed in Table 6. Note that the emerging temperature of 
3.0 keV is consistent with that obtained in spectral fitting 
of the source 1 region in §4.2.1. 

In the first attempt, the source (4C+67.13 nucleus) 
spectrum was fitted by a cosmologically redshifted 
power-law model with the Galactic absorption. The 
background model was fixed and the obtained best-fit 
parameters are listed in Tabled. In this case, the pho¬ 
ton index of the power-law model obtained is quite large 
(r = 2.8), and this model does not provide a very good 
fit (C-stat= 12.4/10) due to a prominent excess around 
1 keV in the data. We have therefore added an addi¬ 
tional thermal component to the fitted source spectrum 
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Source 1 



Source 2-4 



Fig. 10.— Top panels show the observed X-ray spectra and spectral fits (with the abundance fixed at 0.5; red line) for the source 1 (left) 
and the source 2+3+4 (right) regions (see §4.2.1). Bottom panels show the differences between the observed data and the model. 


TABLE 5 

Best-fit model X-ray parameters for the sources 1-4 


Parameter 



Value 



source I 

source 2 

source 3 

source 4 

source 2+3+4 


West 

East 

North 

South 

North & East & South 

Abundance 

0.5 (fix) 

0.5 (fix) 

0.5 (fix) 

0.5 (fix) 

0.5 (fix) 

kT [keV] 

r\ f^o+0.57 
^•‘’^-0.31 

2 04"*"*^'^® 
22 

J..»o_o 26 

i.yo_o 46 


norm'' [I0“^] 

-0.18 

I 

^•'^^-0.16 

i.OZ_o 16 

1.07l“;i® 

^•^'^-0.09 

density [I0“^ cm“^] 

1 71 TO-06 
' ^-0.06 

1 3U+0-07 
i.C)0_o 07 

1 qq + 0.07 
-'-•'^'^-0.07 

l.OSini 

-'■•^'^-0.04 

reduced C-stat/dof 

1.33/37 

1.15/38 

0.79/39 

0.99/34 

1.19/100 


^normalization of APEC model 

(consistently with the analysis results for LINER sources 
in general; see Younes et ah 2011), obtaining this time 
an acceptable fit with C-stat= 4.2/8 (see Table 6) and 
a more reasonable photon index of T = 2.1. The result¬ 
ing temperature of the additional thermal component is 
0.96 keV and the normalization of the power-law compo¬ 
nent is 1.8 X 10“^ photon keV“^ cm“^ s“^. The total nu¬ 
clear 2-10 keV flux in the source rest frame is F 2-10 keV — 
3.4 X 10“^^ ergcm“^ s“^, corresponding to the isotropic 
luminosity of I/ 2 -iokeV — 6.2 x 10^^ergs“^. This is con¬ 
sistent with the bolometric accretion-related luminosity 
Lnuc estimated in § 3 above, assuming the bolometric cor¬ 
rection factor of the order of Tnuc/^2-iokeV ^ 100, as in 
fact expected (Vasudevan & Fabian 2007; Ho 2008). 

5. DISCUSSION AND CONCLUSIONS 

The Abell 578 system provides an interesting insight 
into the widely debated cluster formation processes. Ve¬ 
locity measurements for the member galaxies presented 
by Gomez et al. (1997) indicated that this system is not 
relaxed, and consists of the two interacting sub-clusters. 
Our analysis of newly obtained Chandra data allowed 
us to confirm their results based on ROSAT PSPC data 
for the dominant sub-cluster and to quantify them ro¬ 
bustly. In particular, we found the ellipticity parameter 
for the cluster X-ray atmosphere e 0.3 (i.e., elongated 
east-west), and an offset of about 60kpc between 

the BCG which hosts the 4C+67.13 radio galaxy and 


the centroid of the cluster X-ray emission (inferred from 
the two-dimensional fitting of the surface brightness dis¬ 
tribution of the diffuse emission component). Similar 
displacements have been found in a substantial number 
of galaxy clusters (e.g., Patel et al. 2006; Haarsma et 
al. 2010; Hashimoto et al. 2014). Recently, Lauer et al. 
(2014) demonstrated that the distribution of the pro¬ 
jected spatial offsets is a steep power-law function of 
cluster radius with a median value of about lOkpc. The 
displacement we found for Abell 578/4C+67.13 is there¬ 
fore larger than in most clusters although not as extreme 
as in some other cases (such as, e.g.. Coma cluster with 
rx — 260 kpc), and is also smaller than the Abell 578 
cluster core radius rc — 170 kpc (evaluated here for the 
assumed /3-profile with best-fit index /3 0.7). 

The position of the BCG with respect to the X-ray 
centroid of the cluster is considered a measure of the evo¬ 
lutionary state of a system. In particular, these offsets 
imply deviations from the dynamical equilibrium at ear¬ 
lier stages of a cluster formation via hierarchical merging 
(see, e.g., Martel et al. 2014). In addition, the presence 
of the offsets generally supports the idea that BCGs are 
formed via galaxy mergers in one of the infalling groups 
prior to the cluster assembly (see the discussion in, e.g., 
Merritt 1985; Dubinski 1998; Katayama et al. 2003). In 
this context, we note that the BCG studied here is a 
member of the galaxy pair CGPG 0719.8+6704 (R-band 
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Fig. 11.— Contour plots of temperature versus normalization corresponding to the spectral fits with different abundance values: 0.0 
(top left), 0.5 (top right), 1.0 (bottom left) and free (bottom right). In each panel, the solid and dashed curves denote the 68% and 90% 
confidence levels, respectively, for the northern (black), eastern (red), southern (green), and western (blue) parts of the Abell 578 cluster as 
discussed in §4.2.1 as shown in Figure 9 with the same color-coding. The contours for the combined east+north+south region are indicated 
in magenta. 


magnitudes of 13.6 and 14.8), with the projected sepa¬ 
ration of about lOkpc, and with no or only little young 
stellar populations. The performed optical spectroscopy 
allowed us also to estimate the black hole masses (from 
the observed velocity dispersion) for the two members 
of the pair as Mbh ^ 10^ M© and 3 x 10^ M©, respec¬ 
tively. The analysis of the emission lines indicated more¬ 
over that both galactic nuclei are active, albeit at a low 
level (LINER-type AGN), with the corresponding nu¬ 
clear luminosities of Lnuc ^ 10^^ergs“^. However, only 
the brighter member of the pair is radio-loud. 

The accretion rate for the 4C+67.13 host is rela¬ 
tively low, namely Mace = ^ 0-01 M©yr“^ ^ 

5 X 10“^MEdd5 assuming the radiative efficiency of the 
accretion disk r]d — 0.02 (following Sharma et al. 2007). 
This is comparable to the accretion rate derived by 
Stawarz et al. (2014) for the analogous radio galaxy 
PKSB1358—113 located at the center of the Abell 1836 
cluster. Despite such a limited accretion rate, both AGN 
are able to launch luminous jets surrounded by FRII- 
type radio lobes. Interestingly, the maximum jet kinetic 
power in the source studied here, for the given accretion 
rate as estimated above, reads as Tj/^max) — 3MaccC^ ^ 
2 X 10^^ ergs“^ (see McKinney et al. 2012), which indeed 


is in the range of the jet kinetic luminosities typically de¬ 
rived for “classical” FRIIs (e.g., Machalski et al. 2007). 

We note, on the other hand, that it is quite plausible 
that the nuclear accretion rate in 4G +67.13 has dropped 
significantly since the onset of the jet activity in the 
source, i.e. during the last ^ 1 — lOMyr. And in fact, 
the lack of prominent hotspots at the edges of the radio 
lobes could be considered as an indication for a recently 
quenched/suppressed jet production in the system (due 
to a sudden decrease in the accretion rate), since jet ter¬ 
mination shocks in FRII-type radio galaxies are expected 
to die out almost immediately (i.e., on the timescales of 
^ 10^ — 10^ yr) after the supply of the jet momentum 
ceases (see in this context the discussion in Garilli & 
Barthel 1996, and references therein). If correct, this 
could suggest in general a highly modulated jet duty cy¬ 
cle of BGGs, consisting of short but intense (“FR Il-like”) 
episodes of the enhanced jet activity, separated by the 
extended (“FRI-like”) periods of a quiescence/low-level 
activity. 

In this context, we note that the existing scaling rela¬ 
tions between the jet kinetic power and the lobes’ radio 
luminosity for radio galaxies located at the centers of 
galaxy clusters have been derived assuming a relatively 
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TABLE 6 

Best-fit model parameters for the 4C+67.13 nucleus 


Parameter Value 

Source region 


Model 

power-law 

power-law+apec 

photon index 

2 0 + 0.3 

9 1 + 0.5 

norm^ [ 10 - 6 ] 

4 5“^®'® 
^•^- 0.6 


kT [keV] 

— 

u.yo_0 12 

norm^ [ 10 - 6 ] 

— 

2 q+0.8 

reduced C-stat/dof 

1.24/10 

0.52/8 


Background 

region 

Abundance 

1.0 (fix) 

kT [keV] 

3 

'^•'^-^-0.94 

norm^’ [10“®] 

51 0+9-6 

reduced C-stat/dof 

1.26/6 


^Normalization of the power-law model in units of 
photons keV“^ cm“^ s“^ at 1 keV in the source rest frame. 

^Normalization of the APEC model. 

slow (sonic) expansion of the jet cocoons in a pressure 
equilibrium with the surrounding cluster gas (Birzan et 
al. 2008; Cavagnolo et ah 2010; O’Sullivan et ah 2011). 
When applied to 4C +67.13, these relations give in par¬ 
ticular I/j/(rad) ^ (2 — 7) X lO^^ergs”^. The main as¬ 
sumption involved here may however be justified only 
in the case of FRI systems, but not FRII-type radio 
sources which likely expand with supersonic velocities; 
with such a rapid expansion, leading to the formation 
of bow-shocks in the surrounding intracluster medium, 
the jet kinetic power should be therefore larger than the 
corresponding value of Tj/(rad)- 

Our Chandra data do not provide any direct evi¬ 
dence for the presence of a shock driven in the center 
of Abell 578 by a supersonic expansion of the 4C +67.13 
lobes, but this may only be due to the limited photon 
statistics. Instead, our detailed analysis of the available 
X-ray maps did reveal that the cluster gas in the vicin¬ 
ity of 4C+67.13 seems to be compressed by a factor of 
about 1.4 (i.e., the gas number density increased by this 
factor) and heated (from kT 2.0keV up to 2.7keV), 
consistent with a weak, Aish ^ 1-3 shock driven by the 
expanding jet cocoon. Analogous weak shocks have been 
found in much deeper Chandra and XMM-Newton ex¬ 
posures of several rich clusters (e.g., Nulsen et al. 2005; 
Fabian et al. 2006; Wilson et al. 2006; Gitti et al. 2007; 
Wise et al. 2007; Simionescu et al. 2007; Reynolds et al. 
2008; Lai et al. 2010; Blanton et al. 2011; Cavagnolo et 
al. 2011; Croston et al. 2011), and further, often tenta¬ 
tive detections have been reported for a number of other 


(poorer) systems, including PKSB1358—113 (Stawarz et 
al. 2014, and references therein). 

One should keep in mind, at the same time, that the 
aforementioned limited photon statistics following from 
a rather short Chandra exposure of a relatively low- 
luminosity system, precludes us from making definite 
statements on the precise gas density and temperature 
structure in the central parts of Abell 578. In particu¬ 
lar, it is plausible that the elongated X-ray core of the 
cluster, together with the large offset between the posi¬ 
tion of the BCG and the X-ray centroid inferred from 
the modeling of the diffuse emission component, are re¬ 
lated instead to the sloshing of the intracluster medium 
triggered by the ongoing merger processes (e.g., Tittley 
& Henriksen 2005; ZuHone et al. 2010). Only very deep 
follow-up observations of the studied system can resolve 
this issue by revealing directly the surface brightness and 
temperature discontinuities consistent with “cold fronts” 
— expected in the case of a gas sloshing — rather than 
weak shocks (see in this context Markevitch & Vikhlinin 
2007). 

Altogether, the multifrequency data gathered here for 
the Abell 578 cluster reveal an insightful snap-shot of a 
complex system in formation. At its present stage, it ap¬ 
pears that the cluster galaxies and intracluster gas are 
falling onto the center of the cluster potential, evolving 
slowly toward a dynamical equilibrium, while the cluster 
atmosphere is heated gently by the radio jets/lobes pro¬ 
duced with high efficiency in the low-accretion rate nu¬ 
cleus of the brightest cluster galaxy. The overall cluster 
luminosity To. 5 - 7.0 keV — 2 x lO^^ergs”^ and tempera¬ 
ture kT > 2keV are on the other hand unexceptional, 
being consistent with the X-ray luminosity-temperature 
correlation established for clusters and groups of galaxies 
(see, e.g.. Sun 2012). 
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